Self-induced decoherence approach: Strong limitations on its validity in a simple spin 

bath model and on its general physical relevance 
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The "self-induced decoherence" (SID) approach suggests that (1) the expectation value of any 
observable becomes diagonal in the eigenstates of the total Hamiltonian for systems endowed with 
a continuous energy spectrum, and that (2) this process can be interpreted as decoherence. We 
evaluate the first claim in the context of a simple spin bath model. We find that even for large 
environments, corresponding to an approximately continuous energy spectrum, diagonalization of 
the expectation value of random observables does in general not occur. We explain this result 
and conjecture that SID is likely to fail also in other systems composed of discrete subsystems. 
Regarding the second claim, we emphasize that SID does not describe a physically meaningful 
decoherence process for individual measurements, but only involves destructive interference that 
occurs collectively within an ensemble of presupposed "values" of measurements. This leads us to 
question the relevance of SID for treating observed decoherence effects. 

PACS numbers: 03.65.Ta, 03.65.Yz 



INTRODUCTION 



In a series of papers [ll H, B IJ, IS |a 13 IS IE llfl llll 111 
IT^ IT^ . the authors claim to present a "new approach 
to decoherence" 0, termed "self- induced decoherence" 
(SID). Their main assertion is that, for systems endowed 
with a continuous energy spectrum, the expectation value 
of an observable will become diagonal in the eigenbasis of 
the Hamiltonian of the system, and that this effect can 
be viewed as decoherence. 

The basic idea underlying SID goes back to well-known 
arguments in the context of quantum measurement and 
the theory of irreversible processes [islll^ ll^ IT^ IT^ l20l 

I21I . It rests on the observation that a superposition 
of a large number of terms with random phases in the 
expression for the expectation value of a typical observ- 
able, or for the matrix elements of the density operator, 
leads to destructive interference. The phase differences 
are either due to a random-phase assumption 15], or, as 
in SID, are created dynamically through the time evolu- 
tion factor e*^*/^ associated with each energy eigenstate 
in the superposition. These destructive interference ef- 
fects are then responsible for the diagonalization of the 
expectation value in the energy eigenbasis as described 
by SID. 

However, this process differs strongly from the mecha- 
nism of environment-induced decoherence (EID) [23Ll2^ 
[25, 26, 27, 28, 29, 30, 31]. EID understands decoher- 
ence as the practically irreversible dislocalization of local 
phase relations between environment-selected preferred 
basis states due to entanglement with an environment. 
The approximate diagonality of the expectation value of 
local observables expressed in the preferred basis is only 
a formal phenomenological consequence of the relative 



states of the environment becoming rapidly orthogonal 
during the decoherence process. The fact that SID does 
not require an explicit environment interacting with the 
system motivated the term "self-induced" and was sug- 
gested to circumvent the question of a proper inter- 
pretation of the concept of "observational ignorance of 
the environment" in EID [Hill, SHI. 

This paper pursues two main goals. First, after formal- 
izing the basic idea of SID (Sec.^, we shall discuss the 
question to what extent SID can claim to describe a phys- 
ically relevant decoherence process (Sec. lIII|) . In particu- 
lar, we will argue that, contrary to the claim of its propo- 
nents 0, SID does not constitute a "new viewpoint" on 
decoherence in the usual definition of EID. Second, we 
shall study whether diagonalization of the expectation 
value of random observables in the energy eigenbasis is 
obtained in the context of an explicit spin bath model 
(Sec. lIVp . Dehberately, we have chosen a discrete model 
to investigate the required degree of "quasicontinuity" 
for SID to work as claimed. To anticipate, we find that 
even for bath sizes large compared to what is typically 
considered in EID, no general decay of off-diagonal terms 
is found, unless both the observable and the initial state 
of the bath are appropriately restricted. We explain and 
discuss this result in Sec. and present our conclusions 
in Sec. EH 



II. SELF-INDUCED DECOHERENCE 

The basic formalism of SID as developed in Refs. 0,0, 
1100000113,111111111111 considers an arbitrary 
observable 



6= dE dE' 0{E,E')\E){E'\, 
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(1) 



expanded in the eigenstates \E) of the Hamiltonian H = 
J dE E\E) {E\ with continuous spectrum. In the general 
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treatment, only observables with 

d{E,E') = 0{E)5{E-E') + 0{E,E') (2) 

are considered, where 0{E) and the 0{E,E') are as- 
sumed to be regular functions. The time evolution of the 
expectation value (0)*(t) of O in the pure state = 
e-'-f^*|*o) = / dEe-'^*{E\'^o)\E) (setting ^ = 1) is then 
given by 

(0)*(t) = J dEO{E)\{E\^o)\' 

+ JdE J dE' e-'^^'0{E,E'){E\mo){^o\E'), (3) 

where Ai? = E — E'. For large t, the phase factor e"*^-^* 
fluctuates rapidly with AE, which leads to destructive in- 
terference in the double integral if the multiplying func- 
tion f{E,E') = 0{E,E'){E\^o){'i'o\E') varies compara- 
bly slowly. To formalize this arg ument, SID employs the 
Riemann-Lebesgue theorem '33' , which prescribes that 



(4) 



lim / dzg{z)e''' =0, 



if g{z) is a regular function and integrable (i.e., 
Jdz\g{z)\ < oo). Provided these conditions are satis- 
fied by f{E, E'), it is concluded that 



(t) 



dEOiE)\{E\^o) 



(5) 



for large t. Thus, the off-diagonal terms E ^ E' have 
collectively disappeared, which in SID is interpreted as 
"decoherence in the expectation value." Formally, the 
SID program introduces a "diagonal-equivalent" density 
matrix Orf, 



Pd = J dE\{E\^o)\^\E){E\, 



(6) 



which satisfies (O)p^ = limt^oo (O)p(t). Note that pd is 
only a formal equivalent and is not obtained through any 
dynamical process. Also, expectation values of a nonex- 
haustive set of observables [see Eq. 0] do not uniquely 
determine the density matrix. Therefore, one must not 
derive any conclusions about the possibility for certain 
states of the system from pd- 

To summarize, the main result Eq. (jSJ has been ob- 
tained from two key assumptions: (1) The energy spec- 
trum of the system is continuous; and (2) the coefficients 
used in expanding the initial state and the observable in 
the energy eigenbasis form regular (and integrable) func- 
tions of the energy variable. 

The first requirement of a continuous energy spectrum 
can be viewed as an implicit inclusion of an internal "en- 
vironment" with an infinite number of degrees of free- 
dom. However, any realistic physical system is of finite 
size, and therefore the energy spacing will be discrete. An 



approximate suppression of off-diagonal terms as given 
by Eq. Q should therefore occur also for quasicontin- 
uous energy spectra, i.e., for small but discrete energy 
spacings. 

The regularity assumption (2) is crucial, since it en- 
sures that the phase factors e'^^* are able to lead to the 
required destructive interference of the expansion coeffi- 
cients for large times. However, especially in the realistic 
case of systems of finite size where the expansion coeffi- 
cients will be a finite set of discrete values, this condition 
will not hold. It is therefore important to understand the 
physical meaning and the consequences of a violation of 
this assumption. 

Note also that the strict mathematical limit i ^ oo 
employed in the Riemann-Lebesgue theorem, Eq. |0J, is 
not physically meaningful, and approximate suppression 
must therefore occur already over finite time scales, as in- 
dicated in Eq. (|SJ). Also, for the realistic case of only qua- 
sicontinuous (i.e., essentially discrete) energy spectra, no 
conclusions about an "irreversibility" of the decay should 
be derived from the limit t — > oo (as it is done, for exam- 
ple, in Ref. [2I p. 88]), since the off-diagonal terms will 
return to their initial values within a finite recurrence 
time scale. 

The issues outlined above will be illustrated and in- 
vestigated in the context of a particular model system in 
SecHvl 



III. DOES SID DESCRIBE DECOHERENCE? 

Despite the fact that SID and EID share the term "de- 
coherence" in their name, we shall demonstrate in this 
section that their foundations, scope, and physical im- 
plications are fundamentally different.^ Keeping these 
differences in mind is very important for a proper inter- 
pretation of the study of the bath model described in the 
following Sec. IIVI 

As already briefly outlined in the Introduction, the 
standard a ppr oach of environmental decoherence j23l |2J, 

m m mil m us im describes the consequences of 

the ubiquitous interaction of any system with its envi- 
ronment. This leads to entanglement between the sys- 
tem and the environment and singles out a preferred ba- 
sis of the system that is dynamically determined by the 
Hamiltonian governing the interaction. The relative en- 
vironmental states associated with these preferred states 
rapidly approach orthogonality (i.e., macroscopic dis- 
tinguishability) . Phase relations between the preferred 
states that were initially associated with the system 
alone are now "dislocalized" into the system-environment 
combination due to the entanglement, which consti- 
tutes the decoherence process. In this sense, interfcr- 



^ The author is indebted to H.-D. Zeh and E. Joos for drawing 
strong attention to this point. 



3 



ence between the preferred states becomes locally sup- 
pressed, i.e., decoherence leads locally to a transition 
from a superposition to an apparent ("improper" 1^3) 
ensemble. This can be used to define dynamically in- 
dependent relative local wave-function components that 
can be related to local quasiclassical properties, thereby 
mimicking an apparent "collapse" of the wave function 
[23, 26, 27, 28, 29, 30, 31, 34, 31113. 

The interaction between the system and its environ- 
ment, often referred to as a "continuous measurement by 
the environment," is observer independent and can be 
formulated entirely in terms of wave functions, without 
reference to presumed (classical) concepts such as "values 
of observables" and expectation values (see, for example. 
Chap. 2 of Ref . l23l ) . As it has been emphasized fre- 
quently I2J |2^ l32| . the formalism of local ( "reduced" ) 
density matrices and expectation values presupposes the 
probabilistic interpretation of the wave function and ulti- 
mately relies on the occurence of a "collapse" of the wave 
function at some stage (or on the description of an obser- 
vationally equivalent "branching" process in a relative- 
state framework ^6, 27, 28, 29, 30, 31, 34, 35, 36J). 
The approximate diagonalization of the reduced density 
matrix ps = Tr^pgg (describing the probability distri- 
bution of outcomes of measurements on the "system S 
of interest" immersed into an environment £) in the 
environment-selected basis should therefore be consid- 
ered only as a phenomenological consequence of EID, but 
not as its essence (see also Ref. [23, p. 1800]). Given 
an ensemble of results of measurements of a local ob- 
servable Os, the suppression of off-diagonal terms in ps 
can then be related to the approximate diagonality of 
the expectation value of Os in the preferred basis, since 

{Os)psf: = T^^siPsOs). 

In contrast with EID, SID focuses solely on the deriva- 
tion of a suppression of off-diagonal terms (in the energy 
eigenbasis only) in the expectation value of observables 
pertaining to a single undivided closed system; entan- 
glement through interactions between subsystems plays 
no role in SID. As indicated earlier, the damping effect is 
due to destructive interference between a large number of 
terms with dynamically induced phase differences. Thus 
it is only the averaging process contained in the concept 
of expectation values that leads to a disappearance of in- 
terference terms. Individually, each term remains present 
at all times and is not suppressed independently of the 
other terms. The fact that collectively the off-diagonal 
terms may lead to a mutual canceling-out must not be 
misinterpreted as implying that the measurement "out- 
comes" corresponding to these terms do not occur. Thus 
SID cannot pertain to the relevant problem of a loss of 
interference in individual measurements. In view of this 
argument, the concept of the "diagonal-equivalent den- 
sity matrix" pd, as introduced by the SID program [see 
Eq. ©], is rather misleading, since it gives the incorrect 
impression of an absence of interference terms \E){E'\, 
while the corresponding terms in the expression for the 
expectation value are individually present at all times. 



Derivations of a "classical limit" based on pd |3 appear 
to have overlooked this issue. 

While SID rests on the concept of expectation values, 
i.e., of weighted averages over an ensemble of measure- 
ment outcomes, it does not explain the physical origin 
of these outcomes and their ensembles. In contrast with 
EID, SID does not contain a dynamical account of the 
measurement process itself that could motivate explana- 
tions for how measurement outcomes arise (if only, as 
in EID, in an "apparent," relative-state sense). Conse- 
quently, the assumption of an a priori existence of an 
ensemble of measurement outcomes, as it is inherent in 
SID, could be viewed as a particular application of the 
Copenhagen interpretation. One might then argue that 
in this case decoherence would not even be necessary in 
explaining the observed absense of (macrosopic) interfer- 
ence effects. 

Note that EID makes crucial use of the concept of lo- 
cality in deriving a loss of interference, since globally the 
quantum-mechanical superposition remains unchanged, 
as required by the unitarity of the time evolution of the 
total wave function. As frequently emphasized by Zeh 
(e.g., in Refs. [S^, 113) and others (see, for example, 
Ref. this locality" can be grounded in the (nontriv- 

ial) empirical insight that all observers and interactions 
are intrinsically local. On the other hand, the decom- 
position into a "system of interest" and an environment 
that is ignored from an observational point of view, as 
required in EID, and the resulting implication that the 
relevance of environmental decoherence is restricted to 
local subsystems of the total (nonlocal) quantum Uni- 
verse, has been a subject of ongoing critical discussions 
(see, for example, Refs. |23|, ,2£, _32j ) . Furthermore, 
no general rule is available that would indicate where the 
split between system and environment is to be placed, a 
conceptual difficulty admitted also by proponents of EID 
[27L p. 1820]. These issues seem to have motivated the 
attempt of the SID program to derive decoherence for 
closed, undivided systems. 

However, it is important to note that EID has clearly 
demonstrated that the assumption of the existence of 
closed system is unrealistic in essentially all cases [sslls^ . 
Enlarging the system by including parts of its environ- 
ment, as it is implicitly done in SID in order to arrive 
at a quasicontinuous spectrum, will render the closed- 
system assumption even less physically viable: The com- 
bined system will in turn interact with its surroundings, 
and the degree of environmental interaction will increase 
with the number of degrees of freedom in the system. 
Also, since some interaction with the external measuring 
device will be required, the assumption of a closed sys- 
tem simply bypasses the question of how the information 
contained in the ensemble is acquired in the first place. 
Ultimately, the only truly closed "system" is the Universe 
in its entirety, and one can therefore question the physical 
relevance and motivation for a derivation of decoherence 
for subsystems that are presumed to be closed. 

Furthermore, a general measurement in SID would 
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pertain also to the environment implicitly contained in 
the "closed system," posing the question of how this 
could translate into an experimentally realizable situa- 
tion. And even if such a measurement can be carried 
out, its result would usually be of rather little physical 
interest in the typical situation of observing decoherence 
for a particular object due to its largely unobserved en- 
vironment. 

Finally, in SID, suppression of off-diagonal terms al- 
ways occurs in the energy eigenbasis, which can therefore 
be viewed as the universal "preferred basis" in this ap- 
proach. However, this basis will generally not be useful in 
accounting for our observation of different preferred bases 
for the relevant local systems of interest (e.g., s patial lo- 
calization of macroscopic bodies [231 Isfll (40/4?, 'i^ , chi- 
rality eigenstates for molecules such as sugar ^ , 

and energy eigenstates in atoms j4^). Furthermore, the 
energy eigenbasis cannot be used to describe the emer- 
gence of time-dependent, quasiclassical properties. 

In conclusion, not only is the scope of SID more lim- 
ited than that of EID, but the two approaches also rest 
on different foundations. The interpretation of the pro- 
cesses described by these theories is fundamentally dif- 
ferent, even though phenomenological effects of EID can 
manifest themselves in a manner formally similiar to that 
of SID, i.e., as a disappearance of off-diagonal terms in ex- 
pectation values. Any proposed derivations of an "equiv- 
alence" between SID and EID Q, IS can therefore at 
most claim to describe coincidental formal similiarities 
in the context of very particular models, and only if the 
scope of EID is reduced to the influence on expectation 
values. On the basis of our arguments, we question the 
justification for labehng the process referred to by SID 
as "decoherence." 



IV. ANALYSIS OF THE SPIN BATH MODEL 

By studying an explicit model, we shall now directly 
investigate the claim of SID, that terms not diagonal in 
energy in the expectation value of arbitrary observables 
of the system decay if the system is endowed with a con- 
tinuous energy spectrum. We shall also illustrate formal 
and numerical differences in the time evolution of the 
expectation value of local observables that take into ac- 
count only the degrees of freedom of the system S while 
ignoring the environment £ (the situation encountered in 
EID), and global observables that pertain to both S and 
£ (the case treated by SID). However, in view of our ar- 
guments in the preceding Sec. lIIII this should not be mis- 
understood as a side-by-side comparison of SID and EID. 



While expectation values may share formal similiarities 
in both approaches, they also obliterate fundamental dif- 
ferences between SID and EID that lead to very different 
implications of these expectation values for the question 
of decoherence. 

A. The model and its time evolution 

The probably most simple exactly solvable model for 
decoherence was introduced some years ago by Zurek [2^ . 
Here, the system S consists of a spin-1/2 particle (a single 
qubit) with two possible states |0) (representing spin up) 
and |1) (corresponding to spin down), interacting with a 
collection of N environmental qubits (described by the 
states Iti) and via the total Hamiltonian 

Hse = i(|0)(0| - |1)(1|) f;5.(IT.)(T.| - \U){U\) (g) 

Here, the gi are coupling constants, and = (|Ti)(Ti| + 
|ii)(ii|) is the identity operator for the ith environmen- 
tal qubit. The self-Hamiltonians of S and £ are taken 
to be equal to zero. Note that Hgs has a particularly 
simple form, since it contains only terms diagonal in the 
{|0),|1)} and {IT,}, li.)} bases. 

It follows that the eigenstates of Hs£ are product 
states of the form 10a) = |0)|ti)|i2) ■ ■ ■ ITjv), etc. A gen- 
eral state I^Pq) can then be written as a linear combina- 
tion of product eigenstates, 

N 

Ifo) = {a\0)+b\l))(^{a,\U)+f3^\U)). (8) 

i=l 

This state evolves under the action of Hs£ into 

\^{t))^a\Q)\£o{t))+b\l)\£,{t)), (9) 

where 

N 

\£,{t)) = \£i{-t)) = {a^e^^^-'^'lU) + ^3^e-^'^'/'\U)). 
1=1 

(10) 

The density matrix is 

p{t) = \a\'\0)\£oit)){£„mO\ + \b\'\l)\£,it)){£,itm 
+ ab*\0)\£o{t)){£i{t)\{l\+a*b\l)\£,{t)){£o{t)m, (11) 

and its part diagonal in energy (i.e., diagonal in the eigen- 
states \(p\) of Hgg) is 
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Pd{t) = {\a\'\Q){Q\ + mi){l\) [... + |«in/32p • • • |aivnTi)li2) • • • \U){U\ ■ ■ ■ {l2\{h\ + 



samc-dircction pairing ri2) 

+ (a6*|0)(l| + a*6|l)(0|) [. . . + p^alp^a*, ■ ■ ■ aNf3*N\ii)\i2) ■ ■ ■ \U){iN\ ■ ■ ■ {h\{h\ + •••]• 

" v ' 

opposite-direction pairing 



B. Expectation values of local observables 

Focusing, in the spirit of EID, on the system S alone, 
we trace out the degrees of freedom of the spin bath in 
the density operator pse = |*(t))(^'(i)|. This yields the 
reduced density operator 

Ps = TTePs£ = \am{0\ + \b\'\l){l\ 

+ ab*r{t)\0){l\+a*br*{t)\l){0\, (13) 

where the time dependence of the off-diagonal terms 
|0)(1| and |1)(0| is given by the decoherence factor 

JV 

r{t) = m)\£o{t)} - n(l"^l'^^''* + lApe-^'?'*)- (14) 
The expectation value of any local S observable 

N 



logiok(f)l 



J2 

■s,s'=0,l 



(15) 



is then given by 

{Os)^{t} = ^^^^■ssipOs) = Trs{psOs) 

= lapsoo + Ifepsii + 2Re[ab*swr{t)] .(16) 

We can formally rewrite r(t) as a sum, 

r(i) = 5]|(vI/„|</„)|V^^*, (17) 



where the sum runs over all eigenstates |0a) of the total 
Hamiltonian Hgg, with eigenvalues E\. 

A concrete illustration for the time dependence of r(t), 
Eq. (|14|l . for two different bath sizes is shown in Fig. ^ 
We see that \r{t)\ decays quickly by several orders of 
magnitude and then continues to oscillate about a very 
small mean value. Thus, for local observables, terms cor- 
responding to interference between the two S states |0) 
and |1) become quickly and strongly suppressed. 



C. Expectation values of global observables 

An arbitrary global observable O = Oss can be 
written as a linear combination of the form O = 




FIG. 1: Plot of logj^o |r(t)|, with the decoherence factor r{t) 
given by Eq. (I14II . for two different bath sizes A'^ = 20 and 100. 
Fast decay of r{t), corresponding to local decoherence, is ob- 
served, and the degree of decoherence is seen to increase with 
TV. The squared coefficients \oii\^ and the couplings gi were 
drawn from a uniform random distribution over the intervals 
[0, 1] and [— 7r,7r], respectively. 



Oaa'|0a)(</'a'|; where the \(j)x) are product eigen- 
states of the total Hamiltonian HgSi Eq. O. Explicitly, 



O - E(^o()'|0)(0| + 4i |0)(1| -l-4^/|l)(0| + 4^|1)(1 



N 



(4f iT.)(T.i+4T^iT.)(i.i+4T'ii.)(T.i+4T^u.)a.i)- 



{ij) 



Aij) I 



(18) 



Since O must be Hermitian, sqo, su, elV, and e; are 



real numbers, and sqi = (sio)*: = (^u)*- To keep 
the notation simple, we shall omit the sum over j (and 
thus the index j) in the following. 

The expectation value of O in the state |5'(t)}, Eq. |(^, 



,(0 
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1=1 

^ i=i ' 

= (lapsoo + |6|'sii)ro(t) + 2Re[a6*siori(t)]. 



The special case of the expectation value of local ob- 
servables, as considered in the preceding Sec. II V Bl can 
easily be recovered by remembering that tracing out the 
degrees of freedom of £ is equivalent to choosing all coef- 
ficients ej*| — e^^j^ — 1 and e^*|' — = 0, which yields 
ro(t) = 1 and ri(t) = r{t) [see Eq. (CH)], in agreement 
with Eq. (Unj . 

Suppression of terms in (O)^^^-) that are not diagonal in 
the energy eigenbasis would be represented by the vanish- 
ing of all time-dependent terms in the above expression, 
i.e., 

N 

{d)a = (|apsoo + |5psii)n(l«^lM? + l/^^lM?) 

i=l 

+2 Re I ab*sw J]^ 2 Re{a*P^ef^) \ 
^ i=i ^ 

= (|a|2soo + |fe|'sii)r^ + 2Re(a6*sior^), (20) 

because we can easily show that {0)d = Tr(prfO), where 
Pd, Eq- (|12|l . is the part of the density matrix that is 
diagonal in the eigenstates of the total Hamiltonian. We 
also see that {0)d = Tr(pOrf), where 

N 

da = (.oo|0)(0| + .n|l)(l|) (g) (4?|T.)(T,| + 4?li^)(i^l) 

1=1 

N 

+ (5oi|0)(i| + .io|i>(0|) (g) (4t IT.Xi.l + 4t U')(T'I)- 
i=i 

(21) 

is the part of O diagonal in energy. Thus, as expected, 
diagonality of {0)q,(^t) in energy can also be characterized 
by the presence of only those product expansion coeffi- 
cients that are contained in Od- 

The form of the two product terms ro(i) and Ti{t) is 
similar: They only differ in the order of the pairing of 
the product expansion coefficients with the exponential 
factors. Also, since the coefficients Sjji are independent, 
diagonalization in energy will in general require that indi- 
vidually Tjit) -» rjj and Ti{t) Tf for large t. We can 
therefore restrict our following analysis to ro(t) alone. 
(We shall also omit the subscript "0" in the following.) 



(19) 

I 

First of all, let us rewrite T{t) as a sum of 4^ terms, 
r(i)-^c;,e'^-*. (22) 

A 

where the cx represent products of expansion coefficients. 




Here the sets /^(A) specify over which indices i each 
product runs, namely, they are subsets of the set / = 
{l,...,iV} of all integers between 1 and N such that 
Ufe4(A) = / and HklkW = 0- The total energy Ex 
associated with each term in the sum, Eq. H22|l. is 

Ex= 9^- J2 9^■ (24) 

We choose the index A such that i?A-i < Ex < Ex+i for 
all A. Clearly, Ex = whenever /3(A) = /4(A) = [i.e., 
if /i(A) U /2(A) — I], canceling out the time dependence 
of the associated product term in the expression for r(i). 
Thus, we can split T{t) into a time-independent and a 
time-dependent part, 

m - 5]cA + ^CAe»^^*^r'^ + A(t), (25) 

A A 

where now the first sum runs over all A for which /i(A) U 
/2(A) — I, while the second sum runs over all A for which 
/3(A) U /4(A) ^0. 

Diagonality in energy would require A(t) — > as 
t — > 00. Written this way, we see that A{t) is formally 
similiar to the function r{t) derived for local observables, 
Eq. (|17|l . This might not come as a surprise, since also the 
expression for r(t) can be derived from the calculation of 
an expectation value of an observable, namely, that of the 
local observable Or = (|0)(1| + |1)(0|) (gj^^ Ik that mea- 
sures the degree of local interference between the S states 
|0) and |1). However, in the case of r(i), cx = [((/'aI^'o)!^ 
is a product of TV real and non- negative coefficients [ail^ 
and |/3i|^ = 1 — \ai\^, while the cx of Eq. contain 
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cross terms of the form aiP* and q;*/?^, arbitrary real co- 
efficients e|*| and ej^^ , and arbitrary complex coefficients 

We expect this difference to have strong influence on 
the time evolution of A{t) vs that of r{t). The destruc- 
tive interference needed to obtain suppression of the off- 
diagonal part of the expectation value relies on the idea 
that, when a function f(z) is multiplied by a phase factor 
e'^* whose variation with z is much faster than that of 
f{z), neighboring values f{z) and f{z + Sz) will have si- 
miliar magnitude and phases, but will be weighted with 
two strongly different phase factors, which leads to an 
averaging-out effect in the sum 

In our case, writing 

A{t)=Y,rxe'^'e'^'', (26) 

A 

with r\ = \cx\, the phases ip\ will in general vary very 
rapidly with A and, thus, with Ex- This is a consequence 
of the fact that the c\ = rxe^"^^ are composed of products 
of coefficients, such that changing a single term in the 
product will in general result in a drastic change in the 
overall phase associated with the c\. (The variation in 
magnitude among the cx can be expected to be compara- 
bly insignificant for larger N .) Such discontinuous phase 
fluctuations are absent in the formally similiar function 
r{t), Eq. H14|l . since there only the absolute value of the 
coefficients ai and Pi enters. Note that the impact of 
the phase fluctuations cannot be diminished by going to 
larger t, since the 27r periodicity of phases implies that 
the effect of a phase difference between terms A and A -I- 1 
induced by g^{E\+i-Ex)t ^jj^ average be similiar to that 
induced by e^'^'^^+i"'^^-' for all (larger) values of t. 

We anticipate the described phase-variation effect to 
counteract the averaging-out influence of the multiplying 
phase factor e~*^*, and to thus make it more difficult, 
if not entirely impossible, for A(i), Eq. H25|l. to converge 
to zero. On the other hand, if the average difference 
between the phases associated with the individual coef- 
ficients is decreased, we would expect that the rate and 
degree of decay of h.{t) will be improved. 

D. Numerical results for the expectation value of 
random global observables 

To check this prediction and to generally gain more 
insight into the behavior of {0)^(^t): Eq. 119|) . we studied 
numerically the time evolution of A(t), Eq. H2t)l) . normal- 
ized by its initial value at t = 0, for sets of random ob- 
servables O. Diagonalization of (0)>5(t) in energy would 
then be represented by a decay of A(t) from its initial 
value of one. 

Figure|21shows three typical examples for the time evo- 
lution of logj^o ^(0 fo'" ^ fixed bath size of = 100. All 
couplings gi were taken to be random real numbers be- 
tween — TT and TT. To investigate the influence of phase 



fluctuations of the ca, Eq. (|^ . we considered three dif- 
ferent cases for selecting the coefficients a^, /3i, e||\ ej|\ 

and e||\ i.e., for choosing the initial state of the envi- 
ronment and the observable. In the completely random 
case (A), the coefficients a^, /3i, and e^'j' were taken to be 
random complex numbers, with magnitudes and phases 
drawn from a uniform distribution over the intervals [0, 1] 
and [0, 27r], respectively (and such that \(3i\'^ = 1 — la^p). 

Similiarly, the coefficients ej*| and ej^ were random real 
numbers drawn from a uniform distribution over the in- 
terval [—1,1]. In the second case (B), the initial state 
of the environment was prepared such that the phases 
of the ai and Pi were restricted to the interval [0,7r/2]. 
Also, only observables with non-negative values of e"^}^ 

and el^i^ were considered, such that sign reversals of cx 
due to a change of product terms containing these coef- 
ficients were prevented. Finally, in the third case (C), 
only the absolute values of the ai, Pi, ejp and ej-^^ 
were used, which implies that the cx fluctuated only in 
magnitude. 

We observed a drastic influence of the range of phases 
and signs associated with the individual coefficients ai, 
Pi, ej|, e||, and e^j^j, on the evolution of A(i). In the 
special case (C) of all coefficients being real non-negative 
numbers, A(i) exhibited a consistently strong and fast de- 
cay behavior, similiar to the decay of the function r(t), 
Eq. (|14|l , describing suppression of off-diagonal terms for 
local observables (see Fig. In the intermediate case 
(B), with restricted phases and signs, the degree of de- 
cay of A(t) was decreased, while the decay rate stayed 
roughly the same. In the general random case (A), in 
which no restriction on the spread of phases and on the 
signs of the coefficients was imposed, the time evolution 
of A{t) was observed to be sensitive to the particular set 
of random numbers used for the coefficients ai. Pi, e^^, 

ej*j', and e'"^^ in each run. For some of the sets, A{t) was 
seen to lack any decay behavior at all. In other cases, 
the baseline of oscillation was located below zero, indi- 
cating a very weak damping effect, albeit with the peaks 
of the large-amplitude oscillation frequently reaching val- 
ues greater than zero. 

These results show that, for the bath size studied here, 
a consistent occurence of a decay of A(t) hinges on the 
phase restrictions imposed on the coefficients describing 
the observable and the initial state of the environment. 
If these restrictions are given up, the time evolution of 
A{t) and any occurence of a (comparably weak) decay 
will exhibit strong dependence on the particular set of 
values chosen for the coefficients. 

However, it is important to realize that the assumption 
of a restricted initial state of £ is not only unrealistic, 
since the environment is typically uncontrollable, but it 
will also lead to a circular argument when aiming at a 
derivation of a universal decay effect. This is so because 
any restriction would require an appropriate preparation 
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FIG. 3: Time evolution of log^Q A{t) for iV = 100 bath spins 
see Eq. lESJ] when no restrictions on the initial state of the 
environment are imposed. Such restrictions are physically 
unrealistic and require a preparing measurement on the unre- 
stricted environment, which would in turn be in conflict with 
the desired generality of the derivation of decay effects. No 
collective decay of off-diagonal terms is observed, regardless 
of any restrictions imposed on the observable. 
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FIG. 2: Time evolution of log^o A{t) [see Eq. ^] for iV = 100 
bath spins and three different random observables and ini- 
tial states of the environment. The function A{t) quantifies 
the time dependence of the terms in the expectation value 
(0)„(t) [Eq. Gnj] that are not diagonal in energy. Suppres- 
sion of these terms is represented by a decay of A(i) from its 
initial value of one [i.e., logjQ A(f) — > — oo]. It is observed 
that in the general case (A) of completely random observables 
and initial states of the environment, collective decay of off- 
diagonal terms does, in general, not occur. However, if the 
phases of the coefficients describing the observable and the 
environment are moderately restricted (B) or all set equal to 
zero (C), decay of increasing strength is found. 



of the initial state through a measurement on the entire 
£, which implies that suppression of off-diagonal terms 
would then in general be absent for the observable cor- 
responding to this measurement, if the restriction of the 
initial state of £ is relevant to the occurence of the sup- 
pression. Consequently, the and (3i must be allowed to 



possess arbitrary phases. Then, 



(i) 

are always paired with the and f3i in the expression 
for the c\ that make up A{t) [see Eq. (I23|l ]. we anticipate 
that giving up phase restrictions on the and Pi will 
render the restrictions imposed on the O-coefficients less 
effective, if not entirely irrelevant, in bringing about a 
decay of A{t). 

To study this prediction, in Fig. 13 we show a repre- 
sentative plot of A{t) using only the absolute values of 



-It ■ 



but with the £ co- 



and 



efficients and Pi possessing random phases between 
and 27r. We found that decay is either entirely absent 
or strongly diminished in strength, despite the fact that 
the strongest possible restriction on the phases and signs 
of the O coefficients is imposed. Similiar to the case of 
completely random coefficients, the behavior of A{t) was 
observed to depend crucially on the particular set of ran- 
dom numbers chosen for the coefficients. These results 
lead us to conclude that a universal decay of off-diagonal 
terms does not occur for the studied bath size and time 
scale. 

To be sure, SID is based on the assumption of a qua- 
sicontinuous energy spectrum and very long time scales, 
corresponding to "sufficientfy large" N and f (the existing 
derivations ofSID2,I,", 5,"6,7f|lli[l3[il[ll[ll[l|| 

even assume the strict limits iV — > oo and t — > oo, in 
order to allow for a direct application of the Riemann- 
Lebesgue theorem) , while so far we have only considered 
relatively modest values for these parameters. However, 
since we know from Fig. ^ that for expectation values of 
local observables, strong and fast decay of off-diagonal 
terms is obtained for the value of N and over the time 
scale used in the plots shown in Fig. 12 it is clear that, if 
a general global disappearance of interference terms is to 
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FIG. 4; Example for the time evolution of logjQ A{t) [see 
Eq. 1261 1 using a random observable and random initial bath 
state, for bath sizes between A'' = 10^ and 10^ and a long 
time scale t = 0-10®. No connection between the size of the 
spin bath and the occurrence and the degree of damping is 
observed. Therefore, no consistent collective decay of inter- 
ference terms occurs. The increased time scale is seen to be 
irrelevant. 



occur in our model, it will require a much larger number 
of environmental qubits and/or longer time scales than 
typically considered for local observables. 

Accordingly, in Fig. 0] we show a typical example for 
the time evolution of A{t) over the time scale t — 0- 
10^ for the case of a completely random observable and 
initial state of £, using comparably large bath sizes N 
between N = 10^ and 10^. We observed that even for 
these values of N, no consistent occurrence of a decay 
became apparent. In particular, no generally valid direct 
correlation between the value of N and the time evolution 
of A{t) was visible. Instead, it was again the particular 
set of random numbers included in the computation of 
A{t) for a given value of A'' (but not to the size N of 
the set itself) that determined whether the baseline of 
oscillation of A{t) was located above or below the zero 
line. In agreement with analytical predictions in the pre- 
ceding section, we also found that the choice of a longer 
timescale is irrelevant, since neither the baseline nor the 
amplitude of oscillation changed significantly over the in- 
vestigated time interval after a comparably short initial 
period. Furthermore, we observed that even if A{t) "de- 
cayed" for a particular set of random numbers, the func- 
tion sustained a large-amplitude oscillation whose peaks 
often attained values much larger than the initial value 
of A(t). 

Our results show that, in general, for the bath sizes 
and time scales studied, destructive interference of off- 
diagonal terms in the expectation value expressed in the 
energy eigenbasis [as quantified by A{t), see Eq. (|^ ] 
does not occur in our model. Instead, the time evolution 
of A(t) is simply determined by the particular random 
numbers used to describe the observable and the initial 
state of the environment. Therefore, no general suppres- 
sion of interference terms can be inferred. 



V. DISCUSSION 

The process described by SID appears to be neither 
formally nor conceptually nor physically related to the 
decoherence mechanism in the standard sense of envi- 
ronmental decoherence. EID accounts for the absence 
of interference from the perspective of the local (open) 
system by describing interactions with an environment 
in quantum-mechanical terms of wave-function entangle- 
ment. In contrast, SID describes dynamically induced 
destructive interference between time-dependent terms in 
the expression for expectation values. SID does not, how- 
ever, explain the physical origin of the measurement out- 
comes and their probability-weighted ensembles needed 
to define the expectation values. Even if this purely phe- 
nomcnological basis of SID is accepted, the described pro- 
cess has no bearing on a loss of coherence in individual 
measurements, since it is only a consequence of averaging 
over a large number of measurement results. This is in 
fundamental contrast to EID, where each measurement- 
like interaction leads to a dislocalization of interference 
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and thus, locally, to a disappearance of interference. 

The main result of our study of the spin bath model is 
the finding that the destructive interference predicted by 
SID will in general fail to occur in our model even for bath 
sizes and over time scales much larger than typically con- 
sidered in treatments of the same model in environmental 
decoherence. The source of this failure lies in the random 
relative phases associated with the individual initial bath 
spin states and the expansion coefficients of the observ- 
able. The resulting discontinuous phase fluctuations in 
the coefficient function ca, as defined in Eq. H23() . coun- 
teract the supposed averaging-out effect of the dynamical 
phase factors e'^' in a way that is, due to the 27r period- 
icity of the phase, effectively independent of the value of 
t. 

Even when the bath size is increased, the function c\ 
remains a set of discrete values with discontinuously vary- 
ing phases. This can be explained by noting that, while 
the total energy is a sum of the energies of each sub- 
system, such that enlarging the number of contributing 
subsystems will in general lead to an improved quasicon- 
tinuity of the energy spectrum, the 2tt periodicity of the 
phases implies that the degree of phase discontinuity of 
the Ca will not be diminished by increasing the number of 
subsystems. It is therefore unlikely that a consistent de- 
cay behavior could become apparent for spin baths much 
larger than those considered here. 

This indicates that it is not the degree of continuity 
of the energy spectrum that represents the determining 
factor for obtaining destructive interference. Rather, it 
is the discrete nature of the model itself that seems to 
lead to difficulties. Only if restrictions are imposed on 
both the measured observable and the initial state of the 
environment, a consistent and general suppression of off- 
diagonal terms can occur. But, as we have argued, the 
corresponding preparation of the initial state of the envi- 
ronment is physically unrealistic and renders the deriva- 
tion of a universal decay effect circular. 

We conjecture that the diagonalization of the expec- 
tation value, as described by SID, is likely to fail also in 
other systems composed of discrete individual subenti- 
ties. For, in such models, the relevant function will typi- 
cally be represented by a large product of discrete expan- 
sion coefficients, similiar to the ca of our model, whose 
discontinuous phase fluctuations will again be likely to 
counteract the averaging-out influence of the dynamical 
phases. It is therefore clear that the seemingly innocu- 
ous mathematical requirement of regularity and integra- 
bility of the coefficient functions (see Sec. is far from 
"valid in all relevant cases" where the condition of a suf- 
ficiently continuous energy spectrum holds. The sugges- 
tion to approximate such discrete functions by a continu- 
ous function through interpolation 2] does not appear to 
be viable, since the interpolated function would describe 
a physically different situation. 

On a general note, it is also important to realize that 
dynamical phases are correlated. Thus one could always 
construct an observable for which the initial phases of 



the coefficients seem completely random, but are in fact 
chosen such that recurrence of coherence will show up 
within a finite time interval, thus disproving the claimed 
universality of SID without any further argument. 



VI. SUMMARY AND CONCLUSIONS 

We have investigated the two main claims of the "self- 
induced decoherence" approach, namely, (1) that expec- 
tation values of observables pertaining to a closed sys- 
tem become diagonal in the eigenbasis of the Hamilto- 
nian, provided the system is endowed with a continuous 
energy spectrum; and (2) that this process represents a 
new way of describing quantum decoherence, and that it 
leads to results equivalent to the standard approach of 
environment-induced decoherence. 

We have evaluated the first claim in the context of a 
simple spin bath model of finite size by studying, ana- 
lytically and numerically, the time evolution of expecta- 
tion values of random global observables. We have found 
that, in general, collective decay of terms off-diagonal in 
the energy eigenbasis does not occur over the large range 
of bath sizes and time scales considered. This result is 
not due to an insufficient quasicontinuity of the energy 
spectrum, but is rather rooted in the randomness of the 
phases associated with the observable and the initial state 
of the environment. Even in the limit of large bath sizes, 
the discrete functions for which destructive interference is 
to be derived do not approach their sufficiently smoothly 
varying interpolated approximations required for the dy- 
namical phase averaging to have an effect. 

These results represent an example for a simple model 
system that, although endowed with a quasicontinuous 
energy spectrum, fails to exhibit the decay of off-diagonal 
terms that would be expected from an extrapolation of 
SID to discrete models in the limit of comparably large 
sizes of the system. Such an extrapolation should be pos- 
sible if the approach is to have general physical relevance. 
We have also anticipated that the decay effect described 
by SID will likely be absent also in other similiar models 
that are composed of discrete subsystems. 

With respect to the second claim of the SID program, 
we have questioned the suggestion that SID represents 
a "new viewpoint" on the theory of environment- 
induced decoherence, since the two approaches are based 
on conceptually, formally, and physically unrelated mech- 
anisms. In particular, we have pointed out the following 
key differences and objections. 

(i) SID does not describe the suppression of inter- 
ference for individual measurements, since interference 
terms in the expectation value are not damped individu- 
ally 

(ii) SID simply presupposes the existence of an ensem- 
ble of measurement outcomes, without giving an account 
of its origin in terms of a physical description of measure- 
ment. 

(iii) The assumption of closed systems is unrealistic. 
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especially for systems containing the many degrees of 
freedom needed to obtain the required quasicontinuous 
energy spectrum. 

(iv) The physical feasibility and relevance of measure- 
ments pertaining to the total system-environment com- 
bination is doubtful. 

(v) Energy as the universal preferred basis of the global 
closed system can usually not account for the different 
observed preferred bases for the local system of interest. 

Our study leads us to two main conclusions. First, it 
points to the need for more precise, physically motivated 
criteria for the occurrence of the destructive interference 
effect described by SID. Most importantly, however, the 



physical interpretation and relevance of this effect need 
to be explained. We suspect that the SID approach 
may have mistakenly interpreted and labeled an unre- 
lated process as "decoherence." 
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